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> ' Abstract 

Gluon fragmentation to heavy J PC = 2 ^ quarkonia is studied herein. We compute 

these D-wave states' polarized fragmentation functions and find that they are enhanced by 
"xT ■ 

large numerical prefactors. The prospects for detecting the lowest lying x Di charmonium 



state at the Tevatron are discussed. 
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One of the outstanding challenges in QCD is to understand the process whereby 
colored quarks and gluons hadronize into colorless mesons and baryons. Until recently, 
the only basis for parton fragmentation intuition came from simple models and empirical 
observations. However within the past two years, important progress has been made in 
understanding hadronization from first principles [0-§J. It is now possible to calculate 
the fragmentation functions which specify the probability for heavy quarks and gluons to 
hadronize into quarkonium bound states starting from perturbative QCD. These functions 
involve nonperturbative matrix elements whose values must still be extracted from exper- 
iment or the lattice. But their dependence upon the quarkonium longitudinal momentum 
fraction z can be calculated to lowest order in the strong interaction fine structure con- 
stant a s {m Q ) and the velocity v of the heavy constituent quark inside the bound state. In 
principle, higher order corrections associated with these two small expansion parameters 
may be systematically evaluated as well. These developments have allowed a range of 
hadronization issues to be explored in a limited but model independent context. 

The first fragmentation functions to be computed from perturbative QCD described 
the hadronization of gluons and heavy quarks into S-wave quarkonium bound states |l|,|J . 
These 0(v 3 ) functions can be used to predict the direct production of r\ c and J charmo- 
nia as well as 775 and T bottomonia at lepton and hadron colliders. More recently, 0(v 5 ) 
P-wave fragmentation functions have also been calculated In this paper, we extend 

the ideas and methods developed in refs. |l|-|5[] to the D-wave sector. 

As L = 2 fragmentation functions start at 0(t> 7 ), they are generally less important 
than those for L = and L = 1 quarkonia. P-wave contributions to Q — > Xq have 
been found to be quite suppressed compared to S-wave terms H, and D-wave effects 
should be even smaller still. On the other hand, gluon fragmentation to Xq is known to 
be phenomenologically significant at hadron machines where the number of gluons in the 
initial state is large. Indeed, the dominant source of high p± prompt J/^'s at the Tevatron 
is gluon fragmentation to x c 's followed by single photon emission |6|-|8|]. This prompt J/ty 
mechanism beats all others by almost two orders of magnitude. It is therefore interesting 
to examine the rate at which gluon fragmentation to D-wave quarkonia occurs as well. 

We will focus in particular upon the production of the lowest lying charmonium state 
with the quantum numbers n = 1, L = 2 and S = 0. This J PC = 2 ^ meson has not yet 
been observed. Quark model predictions for its mass fall within the range 3.81 GeV < 
M < 3.84 GeV which lies above the DD threshold P^12|. But since its parity is odd, 1 1 D2 
cannot decay to DD, for the two spinless mesons would have to emerge in an even parity 
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L = 2 partial wave in order to conserve angular momentum. Moreover, the n = 1 state lies 
below the DDn and DD thresholds. Therefore, 1 1 D 2 charmonium predominantly decays 
to lower cc levels or to light hadrons. Its width is consequently narrow. 

The wavefunction for a physical X D 2 quarkonium can be decomposed into a series of 
Fock state components: 

\ l D 2 Quarkonium) = \(QQ) {1) ) + 0{v)\(QQ) {8) g) + 0(v 2 )\(QQ) {1 ' 8) gg) + ■■■. (1) 

The superscript labels on the QQ pairs indicate whether the heavy quark and antiquark 
reside within a color singlet or octet combination. The leading color singlet Fock component 
in (P contributes to g — > X D 2 fragmentation at 0(v 7 ) via the Feynman diagrams illustrated 
in fig. |l|. These graphs mediate X D 2 production through scattering processes like gg — > 
gg* — > gg 1 D 2 . Such gluon fragmentation reactions dominate over lower order parton 
fusion processes when the energy qo of the incoming off-shell gluon g* is large, but its 
squared four-momentum s = q 2 is close to the bound state's squared mass M 2 ~ (2m Q ) 2 . 
In this kinematic regime, the total cross section for gg — > gg x D 2 factorizes up to Oi^sjq^) 
corrections which we neglect 

°{99 -» 99^2) =2x a(gg -> gg) [ dzD g ^i D2 (z). (2) 
The integrated fragmentation probability is thus simply given by the cross section ratio 

\i n t\ gjgg ~1 99^2) 1 f°°ds f 1 , M 2 ,^ , 2 

dzD g ^ D2 ( Z )= — — =^} M - 2 J o d Z 6(s-—)^\A(g -+ D 2 g)\ . 

(3) 

Higher X D 2 Fock state components also participate at 0(v 7 ) in g — > 1 D 2 fragmentation. 
For example, the S-wave (QQ)^ 1 ' 8 -* pair inside the \{QQ)^^ gg) Fock state in eqn. ([]]) can 
be formed through the hard process g* — ► ^q 1 ' g as pictured in fig. |^. It subsequently 
converts to a physical 1 D 2 quarkonium through the spin-preserving emission of two soft 
gluons. The short distance formation of the bound state takes place at 0(v 3 ) while 

the long distance gluon emissions each cost an additional power of v in the amplitude (15[ . 
The graph in fig. ^| is therefore formally of the same order in the heavy quark velocity 
expansion as the diagrams in fig. |l|. 



1 Nonrelativistic potential models yield the estimate T ~ 0.34 MeV for the total width of the 
l x £>2 charmonium state [13,14]. 
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Unfortunately, we do not know any rigorous way to determine the nonperturbative 
matrix elements associated with the higher 1 L>2 Fock state components. However, previous 
experience with gluon and heavy quark fragmentation to P-wave quarkonia provides some 
guidance. In g — > Xq fragmentation, the 0(v 5 ) contribution from the color-octet term in 
the Xq wavefunction 

| 3 P ,i,2 Quarkonium) = |(QQ) (1) ) + 0(v)\(QQ)^g) + ■■■ (4) 

contains one less power of the short distance fine structure constant a s (m Q ) than the 
0(v 5 ) color-singlet term and is numerically more important ||. On the other hand, the 
color-singlet term is much larger than its color-octet counterpart in Q — > Xq fragmentation 
where both are 0(a s (m Q ) 2 ) ||. Since the D-wave diagrams in figs. 1 and 2 are also of the 
same order in a s (m Q ), the P-wave examples suggest it is reasonable to assume that the 
latter graph is numerically small. We will therefore simply neglect it in our analysis. 

We proceed to evaluate the two diagrams displayed in fig. [I] by extending the Feynman 
rules for S and P-wave quarkonium processes derived in ref. ]16[] to D-wave bound states. 
Their sum yields the manifestly gauge invariant amplitude 



iA(g*M - fc 2 (p) + 9b (p')) = -8 A / (g ^ a)3 HZWeMeWeW (p)« (g) 



xe^ T q aP ' T (q + pT(q + P ')P. 



Here N c = 3 denotes the number of colors, g s represents the strong interaction coupling 
and i?2 (0) equals the second derivative of the bound state's radial wavefunction evaluated 
at the origin, i Summing over the final gluon's color and polarization, we obtain the 
squared amplitude 

x q a qW £ e, {q)e v (q) * ^ eg (p)e^ (p) * . 

(7) 



16tt 



The 1 Z?2's polarization tensor is related to its nonrelativistic wavefunction by 



d ' 1 l a lpikH{t\p) = J ^4fl(P)^'(°) (6) 



(2tt)3 a ^nv-^j - Y 8?r c Q/3 



where £ a represents the relative four-momentum between the heavy quark and antiquark inside 
the bound state. The polarization tensor's label h ranges over the helicity levels of the J = 2 
meson. 
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The remaining spin sums for the individual helicity levels of the J = 2 quarkonium may be 



evaluated using the covariant expressions given in ref. fl7|. To the order at which we are 
working, we can also set p^ 1 — zq^ + p± and substitute p^p v — > — i(l — z)(zs — M 2 )g flL ' [l]. 
Then after removing the factor —g' 1L '^2£fi(q)£u(q)* from eqn. (^), inserting the squared 
amplitude into eqn. @, and swapping the order of the s and z integrations, we simply 
read off the polarized gluon fragmentation functions 



D 



>1 D (h=0) 



(z,M) 



16a s (M) 2 |^(0)p 
nN r M 7 



1080z -3 - 234(k -2 + 1680z _1 - 450 + 37z - 2z 2 



+ (1080^" 4 - 2880^" 3 + 2760^" 2 - lUOz' 1 + 190 - 10 z) log(l - 



64« s (M) 2 |^(0)p 

TlNr M 7 



(1-z) 



-360z" 3 + 420z~ 2 - UOz' 1 + 5 + z 



D 



{z,M) 



+ (-360^" 4 + 600^" 3 - 300z -2 + Abz' 1 ) log(l - z) 
32a s (M) 2 |^'(0)| 2 



7vN r 



M 7 



-(1-z) 



180z -3 - 210^" 2 + 30^ _1 + 5 + 2z 



+ (I80z~ 4 - 300z" 3 + 120z -2 ) log(l 



6a) 



and their unpolarized sum 



D 



(z,M) 



80a s (M) 2 \R%(0)\ : 



TTiV, 



M 7 



[3z - 2z 2 + 2(1 - z) log(l-z)]. 



(86) 



The z dependence of this last expression is precisely the same as that of the g — > rj, 
fragmentation function found in ref. [0. 
The functions in eqns. 



sa 



and (pq) are evaluated at the renormalization scale fi = M 
which corresponds to the minimum allowed value for the fragmenting gluon's y/s. They 
may be evolved to higher energies using the Altarelli-Parisi equation 



dD 



, 7 _»id 3 , v a B (fi) f 1 dy ,z v 



where 



Pggiv) = 6 



+ + 2/(1 - y) + ^-3(1 - y) 



(9) 



(10) 



M. 1 y) • y ' ,V1 " '" ' 36 
denotes the gluon splitting function for rif active quark flavors. The \h\ =0,1 and 2 
polarized fragmentation functions for the lowest lying n = 1 1 D2 charmonium state as 
well as their unpolarized sum are plotted in figs. 3a and 3b for p = M and p = 10M 
respectively. The results displayed in the figure are based upon the parameter values 
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M = 3.82 GeV, a s (M) = 0.256 and li^O)] 2 = °- 07 Gey7 0- Comparing the low and 
high energy curves, we see that the weights of all the fragmentation functions are shifted 
toward lower values of z as the renormalization scale increases. This behavior is consistent 
with the general effect of Altarelli-Parisi running upon any fragmentation function. 

An approximate estimate for the rate of prompt Di quarkonia production at hadron 
colliders is given by the product of the total cross section for gluon production and the 
initial integrated fragmentation probability 

80ag(M)|i%(0)| 2 
D ^ {M) ~^N C W— 

In the particular case of l 1 ^ charmonium, we find Dg^ii£> 2 (M) ~ 1.0 x 10 -6 . It is 
instructive to compare this D-wave fragmentation probability with the corresponding S- 
wave result 

D„(M,) = °'' M >»" i 'y,5.3> M -' (12) 

L Vc 

where M Ve = 2.98 GeV, a 8 (M Vc ) = 0.282 and \R (0) | 2 = 0.5 GeV 3 . Although the rate 
for g — > 1 1 I?2 is formally suppressed by four powers of v compared to that for g — > r/ c , the 
probability of the former is enhanced by a numerical prefactor of 80 relative to the latter. 
Gluon fragmentation to the D-wave bound state is therefore larger than one might have 
initially anticipated. 

A more precise prediction for 1 L>2 production may be obtained by folding together 
the Altarelli-Parisi evolved fragmentation functions and the gluon cross section da(pp — > 
g + X)/dp± into the combination 

dajpp^W^ +X) f 1 da(pp^g(p ± /z)+X, t ,) ^ 

dp~± " Jo dp~± ( "' { ] 

This transverse momentum distribution is displayed in fig. [| for 1 1 L>2 charmonium pro- 
duction in the pseudorapidity range \r]\ < 0.6 at the Tevatron. We have used the MRSD0 
parton distribution function evaluated at \i = M±/z = \J M 2 + p\j z to generate the 
differential cross section shown in the figure. 

Since the gluon cross section is a steeply falling function of the main support for 



the integral in eqn. (|13|) lies near z = 1. Looking again at the fragmentation functions 
in figs. 3a and 3b, we see that the h = helicity component dominates over the other 
helicity levels for z close to unity. As a result, the h = differential cross section in fig. [| is 



larger than its \h\ = 1 and \h\ = 2 counterparts. Over the transverse momentum interval 
5 GeV < pj_ < 30 GeV, the ratio of these cross sections is approximately given by 

da( h =V da^^ d a (N=2) 

— : : ~ 1.00 : 0.79 : 0.34. (14) 

dp± dp± dp± 

Recall that the helicity levels for an unpolarized J = 2 state would be populated according 
to 

d a (h=o) Wo-ClfcNi) d a (N= 2 ) 

—, : —, : —, -1:2:2. (15) 

dp± dp± dp± 

Gluon fragmentation consequently induces a sizable X D 2 alignment. This polarization can 
be observed in the angular distribution of photons 

— - oc 1 -0.33 cos 2 6 (16) 
a cost' 

which result from the dominant El radiative transition \ X D 2 ~~ * l 1 -^! +7 |I7|1- A mea- 



surement of this angular distribution would provide a test of the g — > X D 2 fragmentation 
picture. 

The integral of da(pp — > 1 1 D 2 dp± over the transverse momentum range p± > M 
where the results of our gluon fragmentation calculation can be trusted yields 0.8 nb. We 
should stress that this integrated Tevatron cross section value represents a conservative 
lower bound. Higher Fock state contributions to g — > 1 1 -D2, charm quark fragmentation and 
parton fusion processes will all enhance the production of 1 1 D2 charmonia. The D-wave 
state should therefore be produced at a nonnegligible rate. 

Detecting J PC = 2 ^ charmonia will not be simple however. One decay mode which 
might be observable is l 1 !^ — > ijj' ' + 7- This Ml radiative transition connects the x D<i initial 
state to the 3 D\ component of the physical ip' [18]. It is suppressed by the small mixing 
angle which accompanies the L = 2 component. Another possible mode which may be 
experimentally feasible to reconstruct is the following: 



1. 


l x D 2 - 


+ 1*1 + 7 


Br ~ 


0.80 


[13, 14] 


2. 


1*1- 




Br ~ 


0.005 


[19] 


3. 


J/tf - 




Br ~ 


0.06 


[20]. 



Approximate branching ratios for the steps in this decay chain are listed on the right. In 
order to overcome the small value for their product, a data sample corresponding to a 
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large integrated luminosity will have to be collected. The 1 D 2 event rate may then be 
high enough to detect in this channel. 

In conclusion, we have investigated gluon fragmentation to heavy l D 2 quarkonia in 
perturbative QCD. We have found that the polarized g — > X D 2 fragmentation functions are 
enhanced by large numerical prefactors and yield J PC = 2 h mesons which are significantly 
aligned. These results provide motivation to search for these D-wave states in the future. 
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Figure Captions 

Fig. 1. Lowest order color-singlet Feynman diagrams which mediate g — > X D 2 fragmen- 
tation. 

Fig. 2. An example of a higher order Fock component contribution to g — > l D 2 fragmen- 
tation. The ^q 1 ' 8 ^ bound state which is formed at short distances emits two soft 
gluons to become a physical X D 2 quarkonium. 

Fig. 3. g — > 1 1 D 2 charmonium fragmentation functions evaluated at (a) = M and (b) 
\i = 10M. The dotted, dot-dashed and dashed curves represent the fragmentation 
functions for the \h\ = 0, 1 and 2 helicity components of the J = 2 bound state, 
while the solid curve illustrates their sum. 

Fig. 4. Transverse momentum differential cross section for 1 1 D 2 charmonium production 
at the Tevatron. The dotted, dot-dashed and dashed curves represent the cross 
sections for the \h\ = 0, 1 and 2 helicity components of the J = 2 bound state, 
while the solid curve illustrates their sum. 
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